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Air  Force  Research  Lab 


Air  Force  Research  Laboratory 


10  Major  R&D  sites  across  US 
40  Locations  around  the  World 
10  Technical  Directorates 

•  Air  Vehicles  (RB) 

•  Propulsion  (RZ) 

* 

•  Aerospace  Systems  Directorate  (RQ) 


•  5,400  Gov’t  Employees 

•  3,800  On-site  Contractors 


AFRI^i 
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HIGHWAY  14 


Facilities 


Bench-level  Labs 


Altitude  Facilities 

•  From  micro-newtons  to 
50,000  lbs  thrust 


High  Thrust  Facilities 

*  19  Liquid  Engine  stands,  up  to 
8,000,000  lbs  thrust 

*  13  Solid  Rocket  Motor  pads,  up  to 
10,000,000  lbs  thrust 
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Hydrocarbon  Boost 


HCB  establishes  advanced,  modern, 
domestic  LRE  Tech  Base 

—  Required  to  replace  Russian  RD-180  on  EEU 

—  1st  reusable  high  performance  U.S.  HC  engin 

—  Establishes  Ox-rich  staged  combustion 
(ORSC)  tech  base  for  U.S. 

—  Help  sustain  ailing  U.S.  rocket  engine  indust 
tech  development  base 

—  HCB  strongly  supports  SMC/LR  American 
Kerosene  Engine  project 


In-House: 


The  WOWs: 


•  Building  subscale  test  facility  to  mitigate 
combustion  devices  risk 

•  Critical  combustion  research  using  219 
funds 

•  Fuel  thermal 
injector  desig 


•  Design,  build,  test  ORSC  LOx/Kerosene  Liquid 
Rocket  Engine  Tech  Demonstrator 

•  250K-lbf  with  high  Throttle  Capability  (SOTA 
is  2:1)  -  Enables  mission  flexibility 

•  100  Life  Cycle  with  50  cycle  overhaul  (SOTA 
is  20)  -  Exceeds  requirement,  provides  margin 

•  ORSC  is  a  higher  performing  engine  resulting  in  a 
smaller  launch  vehicle  or  an  increase  in  delivered 
payload 
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w  What  is  a  Combustion  Instability  (Cl)? 


Damaged  F-1  engine  injector  faceplate 
caused  by  combustion  instability 

“Combustion  instabilities  have 
been  observed  in  almost  every 
engine  development  effort, 
including  even  the  most  recent 
development  programs” 

-  JANNAF  Stability  Panel  Draft 


•  Combustion  instability  is  an 
organized,  oscillatory  motion 
in  a  combustion  chamber 
sustained  bv  combustion. 

•  Irreparable  damage  can  occur 
in  <1s. 

•  Combustion  instability  caused 
a  four  year  delay  in  the 
development  of  the  F-1  engine 
used  in  the  Apollo  program 

-  >  2000  full  scale  tests 

-  >  $400  million  for  propellants 
alone  (at  2010  prices) 

•  Cl  has  been  identified  as  a 
major  risk  factor  in  the  HCB 
demo  and  future  engine 
development. 
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Risk  Reduction 


Approximate  analysis 


Past  experience 


Growing  CFD 


Candidate 
injector 
designs 


Bombed  testing 

/  \ 


Subscale 
Testing 


0) 
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Systematic  improvements 
systematically  reduce  risk 


Goal  is  to  reach  the 
next  plateau 


Capability  to  model 


Full  Scale 
Design 


Desired  wave  decay 
following  bomb 
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Challenges 


•High  pressures 

-  Supercritical  pressure  with  cryogenic  propellants 

-  Challenging  to  obtain  detailed  data 

•Turbulence  and  Combustion 

-  Unsteady  dynamics  requires  LES  or  hydrid  RANS-LES 

-  Detailed  mechanisms  for  chemical  kinetics 

-  Turbulent  combustion  closures 

•Boundary  Conditions 

-  Simulations  must  include  fuel-ox  manifolds 

•  Data  Processing 
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Overview  of  ALREST 

(Advanced  Liquid  Rocket  Engine  Stability  Technology) 


OBJECTIVE 

•  Develop  advanced  physics-based  combustion 
stability  design  tools  to  reduce  the  risk  of 
developing  combustion  instabilities  in  future  Air 
Force  liquid  rocket  engine  development 
programs. 


APPROACH 

•Fully  coordinate  with  other  national  efforts  to 
conduct  data-centric,  multi-fidelity  model 
development. 

AFRL  9 
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Data-Centric  Model  Development 


Anderson  (Purdue) 

•  AFOSR 

•  NASA  CUI P 

•  ALREST 

•  AFRL 

Frederick  (UAH) 

•  NASA  CUI  P 

•  AFRL 

•  ALREST 
Karagozian  (UCLA) 

•  AFOSR 

Leyva,  Talley  (AFRL) 

•  AFOSR 

•  ALREST 
Cavitt  (Orbitec) 

•  AFRL 

•  ALREST 
Santoro  (PA  State) 

•  AFOSR  (core) 

•  NASA  CUI  P 

•  ALREST 

Yu  (Maryland) 

•  NASA  CUI  P 
Zinn  (GATech) 

•  AFOSR 

Nestleroad  Engin'ng 

•  MDA 
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Experiments 
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Standing  Cl 


Driven 


JetSb  o  • 


Full  Scale  (existing  and  HCB) 
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< - ^ 

Acoustics 
HCB  will  be  heavily 
instrumented  to  provide  Cl 
data  AFRL 
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w  Multi-Fidelity  Model  Development 


Flandro  (GTL) 

•  OSD,  AFRL 
Heister  (Purdue) 

•  AFOSR,  NASA 
Merkle  (Purdue) 

•  NASA,  AFRL,  AFOSR, 
ALREST 

Muss  (Sierra) 

•  AFRL 

Palaniswamy  (Metacomp) 

•  AFOSR,  AFRL,  MDA 
ALREST 

Yang  (PA  State) 

•  AFOSR,  AFRL,  MDA 
Bellan  (J  PL) 

•  AFOSR 

•  ALREST 
Kassoy  (U.  Colo.) 

•  AFOSR 

Priem  consultants 

•  ALREST 
Menon  (GATech) 

•  ALREST,  AFOSR 
Munipalli  (HyPerComp) 

•  ALREST 

Sirignano,  Sideris  (UC 
I  rvine) 

•  AFOSR 
Lynch  (PWR) 

•  ALREST 


Models 
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Near  term  spinoffs 
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Experiments 


ALREST-HFM 

\ 


AAA 

I  f  I 

I 


-Full  physics  to  understand  mechanisms 
and  derive  response  functions 


URANS  with  response 
functions  and  detailed 
geometries 


I 


I 

1 


Response 

Functions 


◄  — 


_ I 
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Analytical  Solution 


|e.g„  Linear 


]•« - 


culery 


Generalized 
Instability  Model 


e.g.,  Galerkin 
Series 


Increased  Cost 


* 
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ALREST-HFM  (AHFM) 


•  ALREST  -  High  Fidelity  Modeling  is  a  six  year 
program  to  develop  high  fidelity  design  tools  for 
combustion  stability 

-  Central  strategy  is  to  take  advantage  of  exponentially 
growing  computational  capability  as  our  fastest  growing 
enabling  tool. 

-  Two  independent  3-year  phases 

*  Selection  for  phase  I  does  not  guarantee  selection  for  phase  II 

•Tools  will  be  validated  against  HCB  data  and 
applied  to  follow-on  engine  programs. 
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Industry  standard  Cl  tools  heritage  empirical  inputs 


Combustion  Stability  Design  Tools 


Current 


iC 
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Industry  standard  Cl  tools 


Combustion  Stability  Design  Tools 


End  of  phase  I 
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Combustion  Stability  Design  Tools 


Future  vision 
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Approach 


Georgia  Tech  Georgia  Tech  HyPerComp 


Turbulence 
Combustion 
Interaction  Model 


LESLIE3D 

•  LES  models 

•  Multiphase  flow 

•  HC  Kinetics 


Turbulence 

Modeling 


Georgia  Tech 

Subcritical 
Spray  Model 


Integration  (HyPerComp) 


Parallel  scaling 
Stability  post 
processing 
Solver 

enhancements 


PWR 


Georgia  Tech 


Phase 

1 


ALREST- 
HFM  1.0 


Enhanced  HC 
Kinetics 

Georgia  Tech 


Supercritical 
EOS  Model 

UTRC 


PWRflow 

Submodels 

PWR 


PWR 

Feeds  Engineering  Models 

•  Chug  code 
• FDORC 

•  Other 


Technology  Demonstration:  GPU,  DG 


Phase 

2 


ALREST- 
HFM  2.0 


Enhanced 

Physics 

Models 


10PD-1 30-044 


Source  code  will  be  delivered  and  maintained 
by  Hypercomp  after  the  contract  ends 
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Marriage  of  PWR  PWRflow  RANS  and  Georgia  Tech 

LESLIE3D  Codes 


Discretization 

Solution 

Advancement 

Schemes 

Types  of 

Grid 

Turbulence 

Modeling 

Kinetics 

Models 

Fluid  Properties 

Turbulence- 

Chemistiy 

Interaction 

Modeling 

Discrete 

Multiphase 

Modeling 

Dispersed 

Multiphase 

Modeling 

L 
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Set  model  for 
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3 

D 

Predictor- 

Corrector 

Temporal  2nd- 

Order 

Unstructured 

Cartesian 

and  K-A 

Hydrocarbon 

Models  for 

Gas  Turbines, 
Ramjets,  and 

LRE’s 

Gases  and 

Mixtures  of  Peng- 

Robinson  Ruids 

with  Detailed 

Species 

Properties 

LES 

modeling 
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flamelets, 
eddy 
breakup 
modeling  to 
LEM  (Linear 
Eddy 
Modeling) 

drop  breakup 

with  interface 

refinement 

Model  with 

drop  breakup 

and 

evaporation 

exercised  on 

gas  turbine 
analyses 

P 

Cell-Based 

Diagonalized 

Unstructured 

Goldberg, 

1- and  2-step 

Perfect  Gas, 

Assumed 

Level  Set 

Lagrangian 

W 

Limiting  with 
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Menter,  and 
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Equilibrium  Air, 

pdf  Model 

model  for  drop 

Droplet 

R 
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Mixtures  of 
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breakup 

Model  with 

fl 
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employed  on 

breakup  now 

o 

(large  density 

Implicit, 

Equation 
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some  selected 

being  tested 

w 

differences) 

Riemann  solvers 

Multigrid  and 

Chunk  G-S 

Qobal  Implicit 

RANS 

Models; 

Anisotropic 

k-  8  Model 

models  for  a 

variety  of 
hydrocarbons 

fuels  and 

hydrogen 
including  RP 

Mixtures  of 

Redlich-Kwong 
and  Feng- 

Robinson  Ruids 

Vulcan  code 

problems 

on  gas 

turbine 

problems 
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AHFM  Dev’t  Team 


\  Pratt  &  Whitney 

A  United  Technologies  Company 

Georgia  1  ©©00@§p 

lech  U  [MlQgjDD'DOOQrOffQg 


United 

Technologies  Research  Center 

yPerComp,  Inc. 
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Four  stages  of  development 


AFRI^i 
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ALREST  verification  suite 


Case  No. 

Description  of  Test  Case  used  for  Verification 

VR-1 

Uniform  Flows  (Run  with  all  available  schemes) 

VR-1.1 

3D  Uniform  Flow  in  rotated  uniform  grid 

VR-1. 2 

3D  Uniform  Flow  in  rotated  non-uniform  grid 

VR-1. 3 

Uniform  Flow  in  a  2-domain  uniform  grid 

VR-2 

Simple  Scaling  Study 

VR-2.1 

3D  Temporal  Mixing  Layer  (TML)  with  light  load 

VR-2. 2 

3D  TML  with  normal  load 

VR-3 

Wave  Propagation  Accuracy 

VR-3.1 

Quasi  ID  Gaussian  pressure  pulse  traveling  in  a  duct  of  variable  area 

VR-3. 2 

Above  with  temperature  variation 

VR-4 

Flame  Test  Cases 

VR-4.1 

Laminar  premixed  methane/air  flame  (phi= 1 , p=  1  to  60  atm,  4-step,  8- 
species,  initial  solution  from  GRI) 

VR-4. 2 

Laminar  premixed  H2/Air  flame  (phi=0.7) 

VR-5 

Boundary  Condition  Test  Cases 

VR-5.1 

Pressure  reflection  from  inflow,  non-reflecting  exit  at  outflow 

VR-5. 2 

Above  with  turbulent  inflow 

VR-5. 3 

Above  with  Calorically  (CPG)  vs  Thermally  (TPG)  perfect  gas  models 

VR-6 

Convection  Test  Cases 

VR-6.1 

ID  Tests  of  wave  speed  with  jump  in  species  concentration 

VR-6. 2 

ID  Shock  tube  problem  with  limiters  and  artificial  dissipation 

VR-6. 3 

ID  Gaussian  pulse  with  different  flux  formulae 

VR-6. 4 

2D  convected  vortex 

VR-6. 5 

ID  Gaussial  entropy  wave 

VR-7 

Temporal  Mixing  Layer 

VR-7.1 

3D,  1  species  Euler  CPG  mixing  layer  model 

VR-7. 2 

2D,  2  species  CPG  model 

VR-7. 3 

Shock  Wave  Test  Cases 

VR-7. 4 

ID  Sod  shock  tube  test  case 

VR-7. 5 

2D  Oblique  shock  Mach  5,  25  deg  wedge 

VR-7. 6 

2D  Richtmyer-Meshkov  Instability 

These  are  the  set  of  “automated  test 
cases”  used  to  verify  code  integrity 
was  maintained  during  code  dev’t 


- 
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2D  blunt  body 


Richtmyer-Meshkov  Instability 


«< 


DENSITY  kg  m3 

1.043e+01 

8.237e+00 

6.038e+00 

3.840e+00 

1.642e+00 
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ALREST  Validation  Cases 


% 
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Magnitude  (psi^/Hz) 


Validation  Simulations 


AFRL  supercritical  data  (validation 
case) 


PSU  LOX/H2  validation  (NASA- 
leveraged) 


Experiment 


itjIlLi. Ji Ll 

2000  4000  6000  8000 


3D  LES 


Single  element  combustion  case 


Purdue  longitudinal  validation  case  wr- 

Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimitecl  bility  Study 


l 


23 


AHFM  Validation  Case: 

Case  3:  CVRC  Longitudinal  Instability  (Yu,  et  al) 


W 


Case  Description: 


Status: 


•  Case  3:  Longitudinal  Instability  for  Single 
Injector 

•  Yu  et  al  Completed  for  Anderson  CUIP  Task 
•Continuous  Variable  Resonance  Combustor 

•  Oxidizer  Post  length= 


Complete;  Analyzing  Results 
Code  Version:  AHFM-al  (LESLIE3D-GA) 
Computer  System  &  No.  of  Cores 
Grid  Size: 

Other  Special  Characteristics: 

*  Mixture  of  Calorically  Perfect  Gases 

*  xx  Methane  -  Oxygen  Kinetics 


Relevance  to  AHFM: 


•  Longitudinal  Instability  for  Hydrocarbon 
Combustion  under  Supercritical  Conditions 


Key  Metric  or  Success  Criteria: 

•  Frequency  and  Amplitude  Growth  of 
Fundamental  Instability  and  Higher 
Harmonic/Secondary  Modes 

•  Mode  Shapes  and  Phase 
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Magnitude  (kP^/Hz)  Pressure  [kPa] 


Pressure  signal 
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Frequency  (Hz) 


Experiment 


LESLIE 


36.85  cm 


Georgia 

Tech 


•  Good  prediction  of 
the  peak  to  peak 
fluctuations 

•  Good  prediction  of 
trends 

•  Frequency  and 
amplitude  slightly  off 

•  200  Hz  and  x2 
respectively 

•  Reason  still  under 
investigation 

•  P0=1.55Mpa> 

P exp=1  -4  MPa 
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Parametric  study  with  Axi-LES 


Good  prediction  of  the  stability 
domain: 

—  L<9cm  and  L>16cm:  strong 
reduction  of  acoustics 

—  L>9cm  and  L<16cm:  unstable 
combustor 

Underestimation  of  the  amplitude 


-t — FWD  2008/11/07 

+ . FWD1  2008/11/21 

■  +-.-FWD2  2008/1 1/21 
■+■■■■  FWD3  2008/11/21 
H — BWD  2008/11/18 
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O  Fixed  Ox  Post 


Effect  of  the  injector  length  on  the 
stability 
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Several-element  transverse  validation 
data  will  come  from  two  phase  II  SBIRs 


Purdue  Jspace 

UNIVERSITY 


Oxidizer  Manifold 


Study  Element 
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CAD  drawing  of  engine  and  manifolds 
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Gloyer-Taylor  Labs’  UCDS  suite  of  tools  applied  to  existing  liquid  rocket  engine  data. 
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Summary 


•ALREST 

-  Nationally  coordinated  data-centric  multi-fidelity  model 
development 

-  ALREST-HFM  is  the  high-fidelity  physics-based  platform 

-  Validated  using  relevant  rocket  data 

-  Results  are  input  into  lower-fidelity  engineering  tools 

•Future 

-  More  sophisticated  physics  models 

-  Improved  combustion  diagnostics 

-  Modular  code  and  model  development 

-  Reduced-basis  model  development 
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